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Optimal  Dynamic  Assignment  for  Low  Earth 
Orbit  Satellite  Constellations 


Alexander  Melin,  R.  Scott  Erwin,  and  VijaySckhar  Chellaboina 


Abstract — In  this  paper  we  investigate  autonomous  task 
assignment  for  a  group  of  low-earth  orbit  satellites,  that  don’t 
necessarily  have  prior  knowledge  of  the  targets  of  interest. 
Specifically,  we  consider  the  optimal  assignment  problem  for 
dynamic  weighted  bipartite  graphs.  First,  we  present  necessary 
and  sufficient  conditions  for  the  existence  of  a  perfect  matching 
in  a  given  bipartite  graph.  Next,  we  present  an  algorithm  to 
construct  a  virtual  graph  based  on  the  original  graph  that 
guarantees  the  existence  of  a  perfect  matching.  These  results 
are  then  used  to  solve  the  optimal  assignment  problem  for 
dynamic  weighted  bipartite  graphs.  Finally  we  apply  this 
algorithm  to  a  constellation  of  low-earth  orbit  satellites. 

Index  Terms — Linear  programming,  dynamic  assignment 
problem,  dynamic  graph  topology,  dynamic  task  allocation, 
graph  theory,  satellite  constellations 


Fig.  1 .  Depiction  of  satellites  tracking  ground  based  targets 


I.  Introduction 

In  recent  years  there  has  been  increased  interest  in  the  use 
of  constellations  of  satellites  and  satellite  formations.  Con¬ 
stellations  of  satellites  are  used  when  global  or  nearly-global 
coverage  are  needed.  Some  applications  include  voice  com¬ 
munication,  satellite  radio,  broadband  networking,  remote 
sensing,  and  laser  communications.  Satellite  formations,  on 
the  other  hand,  employ  close  range  formations  and  the 
functionality  of  the  satellites  depends  on  the  configuration 
of  the  formation. 
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A  lot  of  research  has  been  focused  on  the  control  of 
satellite  formations.  This  is,  in  part,  due  to  the  novelty  of 
the  applications,  such  as  the  Terrestrial  Planet  Finder.  For 
applications  such  as  radar  and  imaging,  the  effective  aper¬ 
ture  of  the  satellite  can  be  increased  without  a  corresponding 
increase  in  the  weight  and  size  of  the  satellite  due  to  the 
replacement  of  the  majority  of  the  structural  elements  with 
virtual  structures.  Much  of  the  research  has  concentrated 
on  controlling  the  geometry  of  the  constellation  and  more 
recently  on  dynamic  reconfiguration  [1], 

Controlling  formation  geometry  lends  itself  to  a  simple 
graph  theoretic  interpretation.  In  general,  a  satellite  constel¬ 
lation  or  UAV  formation  is  considered  as  a  directed  graph 
where  the  nodes  of  the  graph  are  the  satellites  or  UAV’s,  the 
arcs  of  the  graph  are  the  communications  channels  and  the 
physical  position  of  the  nodes  corresponds  to  the  physical 
position  of  the  satellites  or  UAV’s  [2],  [3]. 

On  the  other  hand,  satellite  constellations  such  as  Glob- 
alStar  do  not  rely  on  the  relative  formation  to  provide 
their  functionality,  but  the  specific  orbits  but  the  orbits 
the  satellites  are  placed  in.  These  constellations  are  gen¬ 
erally  used  for  global  radio-frequency  communication.  RF 
communication  can  achieve  global  or  near-global  coverage 
with  a  relatively  small  number  of  satellites  due  to  the  wide 
coverage  area  of  an  individual  satellite  and  overlapping 
footprints. 

Laser  communication  and  most  remote  sensing  satellites 
have  a  narrow  beamwidth  so  simultaneous  global  coverage 
would  require  a  prohibitively  large  number  of  satellites.  To 
achieve  global  coverage  with  narrow  beamwidth  satellites, 
one  solution  is  to  place  the  satellites  into  a  global  coverage 
constellation  with  overlapping  footprints  and  assign  the 
satellites  to  look  at  targets  in  their  footprint  in  some  an 
manner. 

Currently,  for  many  remote  sensing  systems,  the  assign¬ 
ment  of  targets  to  a  satellite  assumes  that  the  target  and 
satellite  positions  are  known  ahead  of  time.  This  allows 
a  schedule  to  be  calculated  ahead  of  time  and  passed  to 
the  relevant  satellites.  In  the  case  of  laser  communications 
this  means  that  a  pre-determined  time  window  is  scheduled 
for  a  connection.  A  more  effective  laser  communications 
system  would  autonomously  reconfigure  itself  to  allow  the 
user  bandwidth-on-demand. 

Performing  this  scheduling  on  the  ground  means  that  the 
system  is  not  quickly  adaptable  to  changes,  in  addition,  de¬ 
lays  are  introduced  by  the  need  to  recalculate  and  retransmit 
the  assignments. 


A.  Problem  Assumptions 

This  paper  proposes  a  method  based  on  the  optimal 
assignment  problem  in  linear  programming  that  allows 
a  group  of  satellites  to  autonomously  change  their  task 
assignments  in  a  distributed  process.  In  order  to  make  the 
problem  more  concrete,  we  will  look  at  a  general  group  of 
satellites  and  make  some  assumptions  on  their  properties 
and  abilities. 

Specifically,  we  will  look  at  groups  of  satellites  whose 
tasks  are  to  acquire  and  track  specific  targets  on  the  surface 
of  the  earth  for  purposes  of  communication,  imaging,  etc. 
The  satellites  must  have  the  ability  to  communicate  with 
each  other.  Since  we  are  concerned  with  high  level  task 
assignment,  the  satellites  are  assumed  to  have  a  stable 
attitude  controller,  and  be  able  to  track  a  point  on  the  earth 
given  the  target’s  position. 

We  will  assume  that  each  satellite  has  the  capability 
to  perform  any  task  assigned  to  it  (i.e.  track  any  target 
that  is  in  range).  Because  of  the  overlapping  nature  of  the 
satellites’  abilities  to  perform  the  tasks,  the  objective  of  the 
dynamic  scheduling  problem  is  to  keep  satellites  from  trying 
to  perform  the  same  task  simultaneously,  while  choosing  the 
tasks  that  maximize  the  overall  efficiency  of  the  group. 

A  purely  decentralized  control  architecture  for  this  prob¬ 
lem  would  require  no  communication  between  satellites, 
and  each  satellite  would  require  knowledge  of  only  its  own 
internal  states.  We  will  assume  an  architecture  which  is 
predominantly  decentralized  but  allows  for  minimal  com¬ 
munication  between  the  satellites,  as  follows:  the  satellite 
must  be  able  to  calculate  three  pieces  of  information  and 
must  transmit  and  receive  one  piece  of  this  information. 

i)  Each  satellite  must  know  all  the  potential  targets  and 
their  locations. 

ii)  Each  satellite  must  be  able  to  determine  which  targets 
are  in  range  and  which  are  out  of  range. 

Hi)  Each  satellite  must  be  able  to  calculate  a  performance 
metric  that  describes  the  efficiency  of  the  satellite 
tracking  each  target  in  its  own  range. 

iv)  Each  satellite  must  be  able  to  transmit  this  vector  of 
performance  metrics  to  all  other  satellites. 

The  control  objective  is  to  maximize  the  time  each  satel¬ 
lite  spends  tracking  targets  over  an  interval  of  time,  subject 
to  the  limitations  of  the  satellites’  abilities.  This  problem 
also  has  many  other  similar  applications  such  as  missile 
tracking  and  can  significantly  reduce  operator  workload,  and 
is  especially  valuable  in  time-critical  military  applications 
[4].  Because  of  the  assumptions  on  the  homogeneity  of 
the  satellites’  abilities  and  the  tasks  to  be  assigned,  the 
problem  can  be  described  by  a  graph  theoretic  approach, 
or  an  equivalent  problem  in  linear  programming  known  as 
the  assignment  problem. 

B.  Dynamic  Graphs  and  the  Assignment  Problem 

The  assignment  problem  in  graph  theory  and  linear 
programming  has  been  well  understood  since  1955  when 


Kuhn  used  a  theorem  by  the  Hungarian  mathematician 
Konig  to  create  an  algorithm,  known  as  the  Hungarian 
method,  that  finds  an  optimal  perfect  assignment  [5].  A 
more  precise  definition  of  an  optimal  perfect  assignment 
will  be  given  later.  Since  this  algorithm  yields  the  exact 
solution  to  this  problem,  later  work  in  this  field  is  concerned 
with  improving  the  speed  of  the  algorithm  [6],  [7].  In  graph 
theoretic  notation,  this  problem  is  a  matching  on  an  edge 
weighted  bipartite  graph  where  one  set  of  nodes  corresponds 
to  the  satellites,  the  other  set  of  nodes  corresponds  to  the 
targets,  the  edges  of  the  graph  connect  the  satellites  and 
targets  that  are  in  range,  and  the  weights  on  the  edges 
correspond  to  a  performance  metric  for  a  satellite  to  track 
that  target. 

The  solution  to  the  optimal  assignment  problem  assumes 
that  the  graph  is  static,  however  this  problem  cannot  be 
described  by  a  static  graph.  The  number  of  satellites  and 
targets  are  constantly  changing,  and  the  edges  of  the  graph 
also  change  as  the  satellites  orbit  the  earth  and  targets  move 
in  and  out  of  range  of  satellites.  Finally  the  edge  weights 
are  changing  as  the  satellite’s  attitude  changes  to  track  a 
specific  target. 

Unfortunately  very  little  is  known  about  dynamic  graphs 
[8],  however,  some  research  is  being  done  in  the  area 
of  combinatorics  and  graph  theory  applied  to  spacecraft 
formation  flying  [1],  [9], 

The  problem  of  controlling  large  groups  of  autonomous 
vehicles  has  been  studied  in  many  forms.  This  is  an  ex¬ 
tremely  complex  problem.  Even  finding  a  universal  model 
for  these  distributed  systems  has  not  been  accomplished 

[10] .  One  related  problem  is  the  scheduling  problem.  The 
most  common  uses  for  this  formulation  are  in  flexible 
manufacturing  systems,  computer  networks,  and  processor 
resource  allocation.  In  general  this  problem  is  NP-complete 

[11] ,  Some  simplified  cases  have  been  solved  in  polynomial 
time,  but  due  to  the  complexity  of  the  problem  heuristic 
methods  are  the  preferred  solution  to  this  problem.  The 
scheduling  problem  consists  of  allocating  tasks  or  jobs  to 
resources  so  that  an  optimality  condition  is  satisfied.  The 
most  common  optimality  condition  is  minimizing  the  total 
job  time  [11],  [12], 

Another  approach  to  controlling  groups  of  autonomous 
vehicles  is  cooperative  control.  This  approach  is  concerned 
with  controlling  the  interactions  between  the  elements  of 
the  system  so  that  a  goal  is  accomplished.  The  tasks  stud¬ 
ied  in  coordinated  control  generally  require  very  complex 
interactions,  such  as  RoboFlag  [10],  Generally  a  ’’best” 
solution  to  these  problems  cannot  be  determined.  Most  of 
the  literature  on  cooperative  control  deals  with  languages 
to  program  the  interactions  between  the  robots  so  that 
the  interactions  are  stable  and  heuristic  strategy  algorithms 
can  be  implemented  to  achieve  the  desired  behavior.  Many 
applications  of  cooperative  control  also  use  decentralized 
control  as  part  of  their  control  paradigm. 

The  problem  most  similar  to  this  problem  is  the  design  of 
network  routing  protocols.  The  main  difference  being,  that 


there  is  a  cost  associated  in  reconfiguring  a  satellite  constel¬ 
lation,  while  reconfiguring  a  network  does  not  involve  any 
loss  of  performance.  Reconfiguration  of  network  is  done 
mainly  to  prevent  channel  overloads  and  reconfigure  in  the 
event  of  a  channel  failure. 

One  main  similarity  in  these  approaches  is  that  all  the 
goals  of  the  problem  require  a  high  level  of  interaction 
between  the  separate  elements  of  the  system  for  the  goal 
to  be  accomplished.  In  flexible  manufacturing,  some  jobs 
must  be  completed  before  others.  In  many  instances  a  single 
agent  cannot  complete  the  goal.  In  cooperative  control 
the  formation  and  interaction  between  the  agents  is  also 
necessary  for  the  completion  of  the  task.  In  autonomous 
formation  flying,  for  example,  the  function  of  the  spacecraft 
depends  on  the  formation  itself. 

The  organization  of  the  remainder  of  the  paper  is  as 
follows.  In  Section  II  we  will  introduce  some  notation 
and  cover  the  concepts  of  graph  theory  required  for  the 
remainder  of  the  paper.  We  also  introduce  the  optimal 
assignment  problem  and  discuss  the  solutions  to  it  that  exist 
in  the  literature,  and  prove  the  existence  of  a  solution.  In 
Section  III  we  introduce  the  concepts  of  the  dynamic  graph 
and  discuss  how  applying  optimal  assignment  problem  to  a 
dynamic  graph  affects  the  behavior  of  the  algorithm. 

II.  Mathematical  Preliminaries 

In  this  section,  we  introduce  notation  and  some  key 
results  necessary  for  developing  the  main  results  of  this 
paper.  Let  R.  (resp.,  Z)  denote  the  set  of  real  numbers  (resp., 
positive  integers),  let  Rnxm  (resp.,  Zraxm  denote  the  set  of 
real  (resp.,  positive  integer)  matrices,  and  let  AT  denote  the 
transpose  of  A. 

Let  G  =  (V,  E )  denote  a  bipartite  graph  where  V  = 
X  x  Y  is  the  set  of  vertices  and  E  is  the  set  of  edges.  It 
is  convenient  to  assume  that  the  elements  of  X  and  Y  are 
enumerated  for  {1, . . . ,  n}  and  {1, . . . ,  to}  so  that  X  and 
Y  may  be  assumed  to  belong  to  Zn  and  Zm,  respectively. 
Furthermore,  E  can  be  represented  by  the  connectivity  (or 
adjacency)  matrix  given  by 

1,  if  {i,j)  is  an  edge 
0,  otherwise 

for  all  i  £  {1, . . . ,  n}  and  j  £  {1, . . . ,  to}.  Next,  let  G  = 
(V.  E,  W)  denote  a  weighted  bipartite  graph  where  (V,  E) 
is  a  bipartite  graph  with  W  £  R"x"  denoting  the  weight 
matrix  associated  with  the  edges.  We  assume  that  IL');  >  0, 
if  Eij  =  1  and  Wij  =  oo  if  E^  =  0. 

A  vertex  i  £  X  (resp.,  j  £  Y)  is  incident  to  an  edge  if 
there  exists  j  £Y  (resp.,  i  £  X)  such  that  (i,j)  is  an  edge. 
In  this  case,  (i,j)  are  adjacent.  A  set  of  vertices  U  C  V 
is  a  vertex  cover  of  G  if  the  set  {(i,j)  :  El;l  =  1  ,i  £ 
{1, . . . ,  n},j  £  {1, . . . ,  to}}  C  U  that  is,  if  every  edge  in 
E  is  incident  with  a  vertex  in  U .  A  minimum  vertex  cover  is 
a  vertex  cover  of  G  with  the  minimum  number  of  vertices. 
The  existence  of  such  a  minimum  vertex  cover  is  easy  to 
establish. 


Two  edges  of  a  graph  are  independent  if  they  do  not 
have  any  common  incident  vertices.  A  matching  of  G  is  a 
set  MCE  which  is  a  subgraph  of  independent  edges  in 
a  graph  G  =  (V,E).  Since  G  is  bipartite  it  follows  that 
the  set  of  incident  vertices  U  of  a  matching  will  be  of  the 
form  X  xY  such  that  X  and  Y  have  the  same  number  of 
vertices.  A  maximum  matching  of  G  is  a  matching  of  G 
that  has  the  largest  number  of  edges  (and  hence  the  largest 
number  of  incident  vertices).  A  perfect  matching  of  G  is 
a  (maximum)  matching  of  G  such  that  the  set  of  incident 
vertices  is  X  x  Y  (see  Figure  2). 


Fig.  2.  A  complete  bipartite  graph,  a  matching  on  a  bipartite  graph  and 
the  equivalent  edge  matrix 

By  definition,  if  a  graph  G  has  a  perfect  matching  M 
then  m  =  n,  that  is  X  and  Y  have  the  same  num¬ 
ber  of  vertices.  Furthermore,  Me  =  MTe  =  e  where 
e  =  [  1  1  1  . . .  ]  .  Note  that  these  conditions  are 

only  necessary  but  not  sufficient  for  the  existence  of  a 
perfect  matching.  Hence,  it  follows  that  if  a  bipartite  graph 
G  =  ( V,E )  has  a  perfect  matching  then  Ee  >>  e  and 
ET e  >>  e.  However,  note  that  these  conditions  are  not 
sufficient  (for  example,  consider  the  graph  given  in  Figure 
3).  The  following  result  provides  a  necessary  and  sufficient 
condition  for  the  existence  of  a  perfect  matching.  First, 
however,  we  introduce  the  following  notation.  Let  Afjj  : 
Rm  x  n  — >  R  be  defined  by 

Md{E)  =  eTmDen  (2) 

where  D  £  Rmxra  is  given  by 


Fig.  3.  A  minimum  line  covering  on  an  edge  matrix 

Note  that  if  m  =  n  then  Md{E)  denotes  the  number 
of  one  on  the  diagonal.  Next,  let  umax  :  Rmxra  — >  R  be 
defined  by 

vmax{E)  =  max  JfD(PmEPn)  (4) 
prneurn,pneuri 


where  Wm  C  Rmxm  and  n„  C  Rnxn  denote  the  set  of 
m  x  to  and  n  x  n  permutation  matrices  [?],  respectively. 

Theorem  2.1:  Consider  a  bipartite  graph  G  =  (V,  E) 
where  E  £  Rmxn.  The  size  of  the  maximum  matching 
is  vmax(E).  Furthermore,  if  to  =  n  then  G  has  a  perfect 
matching  if  and  only  if  vmax  =  n. 

Proof:  Let  M  denote  a  maximum  matching  of  G.  In 
this  case  it  is  easy  to  show  that  there  exists  permutations  of 
X  and  Y  such  that  the  edge  matrix  Mn  of  the  permutations 
Xu  and  Yu  is  such  that  Mn,  =  0,  i  f  j  and  the  size  of  the 
matching  Mu  is  given  by  e„Mne„  =  The  result 

now  follows  from  the  fact  that  there  exist  Pm  £  Hm  and 
Pn  £  n„  such  that  Mn  =  PmMPn  and  ND{PmEPn)  = 
tt'm  Apr  tin  • 

Next,  if  to  =  n  it  follows  that  G  has  a  perfect  matching  if 
and  only  if  vmax(E)  =  n  by  noting  that  a  perfect  matching 
is  a  maximum  matching  of  size  n.  ■ 

Remark  2.1:  It  follows  from  Konigs  theorem  [13]  that 
Vmax(E)  is  the  size  of  the  minimum  vertex  cover.  Hence, 
Vmax(E)  denotes  the  minimum  number  of  lines  (drawn 
across  rows  and  columns)  that  are  sufficient  to  cover  all 
the  l’s  in  E  which  implies  that  vmax{E)  is  the  maximum 
number  of  l’s  no  two  of  which  are  in  the  same  line  (row 
or  column)  of  E 

Remark  2.2:  Although  Theorem  2.1  gives  necessary  and 
sufficient  conditions  for  the  existence  of  a  perfect  matching 
it  is  in  general  more  computationally  efficient  to  use  Konigs 
theorem  to  compute  vmax(E). 

Next,  we  consider  the  virtual  perfect  matching  problem 
where  given  a  graph  G  =  (V,  E)  that  does  not  have  a  perfect 
matching  we  provide  a  method  for  constructing  the  virtual 
graph  G  =  (V,E)  such  that  G  has  a  perfect  matching  and 
G  is  a  projection  of  G ,  that  is,  G  is  a  subgraph  of  G  such 
that  when  all  the  vertices  in  V\V  are  eliminated  (along 
with  their  edges)  then  the  resultant  graph  is  G.  Note  that  if 
E  £  Rmxn  then  by  choosing  V  =  X  x  Y  C  rx^rx" 
such  that  all  the  n  additional  (virtual)  vertices  of  X  and  the 
to  additional  virtual  vertices  of  Y  are  connected  with  every 
other  vertex  of  Y  and  X,  respectively,  it  can  be  shown  that 
G  has  a  perfect  matching.  The  following  result  provides 
a  construction  G  that  has  the  minimum  number  of  virtual 
vertices. 

Theorem  2.2:  Let  G  =  (V,  E)  be  a  bipartite  graph.  Then 
there  exists  G  =  ( V ,  E)  such  that  G  is  a  projection  of  G, 
G  has  a  perfect  matching,  and 

^max(P)  =  (TO  T  Tl)  r/max(P)  (5) 

f  ;u  rt  her  mo  re,  there  does  not  exist  G  =  (V ,  E)  such  that 
G  is  a  projection  of  G,  G  has  a  perfect  matching,  and 

Vmax  (■ E)<  Vmax  (E). 

Proof:  Let  M  be  a  maximum  matching  and  assume 
without  loss  of  generality  that  M,;?  =  0 ,  i  f  j.  Now,  it  can 
be  shown  that  by  adding  n  —  vmax(E)  and  m  —  vmax(E ) 
virtual  vertices  to  X  and  Y  the  resultant  G  has  a  perfect 
matching.  The  minimality  of  vmax(E)  given  by  (5)  can  be 
shown  in  a  similar  manner.  ■ 


Remark  2.3:  Note  that  if  G  has  a  perfect  matching,  that 
is,  vmax(E )  =  m  =  n  is  follows  that  vmax(E)  =  vmax(E ) 
and  hence  G  =  G  (since  G  is  a  projection  of  G) 


Fig.  4.  Augmented  edge  matrix  with  the  new  minimum  line  cover 

Finally,  we  consider  the  optimal  assignement  problem. 
Specifically,  let  G  =  ( V ,  E,  W)  be  a  weighted  bipartite 
graph  that  has  at  least  on  perfect  matching.  The  optimal 
assignment  problem  is  the  problem  of  finding  a  perfect 
matching  that  has  the  minimum  weight ,  that  is,  finding  M 
such  that 

M  =  min  eT  (A4  o  M)e  (6) 

MeM 

where  XI  denotes  the  set  of  all  perfect  matchings  of  G. 
More  generally,  let  G  =  (V,  E,  W)  be  a  weighted  bipartite 
graph  such  that  E  £  Rmxn  which  may  not  have  a  perfect 
matching.  In  this  case,  the  optimal  assignment  problem  is 
given  by 

M  =  min  eT  (W  o  M)e  (7) 

MeM 

where  M  denotes  the  set  of  perfect  matchings  of  G,  the 
virtual  graph  of  G,  and  W  denotes  the  weighting  matrix 
of  G  such  that  Wij  =  0  if  i  or  j  correspondes  to  a  virtual 
vertex  and  Wtl  =  W,j  otherwise. 

Note  that  the  construction  of  G  can  be  performed  using 
the  algorithm  described  in  [16]  and  the  optimal  assignment 
problem  given  by  (7)  can  be  solved  by  using  the  Hungarian 
method  [5],  [13],  [16]. 

III.  Decentralized  Decision  Algorithm 

The  Hungarian  Method  was  meant  to  solve  for  a  mini¬ 
mum  weight  perfect  matching  on  a  static  graph,  but  a  more 
relevant  problem  with  regards  to  the  application  of  interest 
in  this  paper  is  finding  a  minimum  weight  perfect  matching 
as  a  graph  is  dynamically  changing.  In  this  section  we  will 
define  the  problem  of  minimizing  the  total  edge  weight  for 
a  dynamic  graph.  We  will  also  discuss  the  conditions  for 
applying  this  algorithm  to  a  dynamic  graph. 

A  dynamic  graph  is  a  graph  where  the  vertices,  vertex 
weights,  edge  matrix,  and  edge  weight  are  changing  with 
time.  Since  we  are  dealing  with  a  matching  on  a  bipartite 
graph,  vertex  weights  are  not  used  and  can  be  assumed 
constant.  We  will  denote  a  dynamic  graph  by  Q(t)  = 
(V(f),£(f),  W(f))  where  8(t)  and  W(t)  are  nxn  matrices. 

In  the  static  case,  the  edge  weight  matrix  is  a  matrix  of 
positive  constants.  In  the  dynamic  case  W  depends  on  time 
and  the  time-varying  edge  matrix.  This  leads  to  the  general 
definition  of  a  time-varying  weight  matrix  for  a  bipartite 
graph  given  by 

,s) 


where  :  K.  — >  X  C  M+  are  a  set  of  functions  that 

describe  the  weight  function  for  a  particular  edge  of  Q  and 
are  not  necessarily  smooth  or  continuous,  and  1  is  a  closed 
bounded  interval. 

The  choice  of  the  weight  function  is  determined  by  the 
problem  objective.  In  the  case  of  LEO  sensing  or  imaging 
constellations,  the  objective  is  to  maximize  the  time  that 
the  satellites  are  tracking  a  target.  This  is  equivalent  to 
minimizing  the  time  a  the  satellites  spend  moving  between 
targets.  For  this  problem  the  optimal  assignment  problem  is 
to  minimize  the  integral  of  the  total  weight  of  the  matching 
over  an  interval  of  time. 

The  total  weight  of  the  matching  is  given  by 

n  n 

w{m)  (9) 

*=  1  3= 1 

This  leads  to  the  continuous  optimal  assignment  problem 
given  by 

M(t)  =  /  min  w(M)  dt  (10) 

Jt !  M 

„ts  n  n 

=  /  ^EE^vW-^v  (n) 

•>ti  i- 1  j~% 

where  A/f(£,  W)  is  the  set  of  admissible  matchings  at  time 
t. 

Equation  9  shows  that  the  behavior  of  the  total  weight 
function  is  dependent  on  the  matching  and  the  edge  weight 
functions.  For  example,  if  the  matching  does  not  change 
on  the  time  interval  and  the  edge  weight  functions  are 
continuously  differentiable,  then  the  total  weight  function 
will  be  continuously  differentiable  on  the  interval.  On  the 
other  hand,  if  the  edge  weight  function  is  constant  on  the 
time  interval,  the  matching  will  only  change  when  the  edge 
matrix  changes.  This  means  that  the  matching  and  edge 
weight  functions  must  be  “well  behaved”  for  the  integral  to 
behave  nicely. 

In  the  following  we  show  that  if  the  total  edge  weight  is 
bounded  on  the  time  interval,  then  the  solution  is  optimal 
and  the  finite  time  version  of  the  problem  approaches  the 
optimal  solution  in  the  limit.  This  is  equivalent  to  a  perfect 
matching  existing  and  ftJ  (f),  i,j  =  1,2 ,n  bounded  on 

In  actual  implementation,  the  algorithm  will  be  applied 
in  a  discrete  manner  with  the  matching  changing  only  at 
discrete  times.  In  the  following  we  show  that  as  the  time 
between  the  matching  updates  decreases,  the  algorithms 
behavior  approaches  the  continuous  time  solution. 

We  will  start  by  defining  a  discrete  time  version  of  the 
optimal  assignment  problem.  We  will  allow  the  matching 
to  be  made  at  discrete  times  in  the  time  interval.  Let  P  = 
{tk}l=1  be  a  sequence  of  times  that  partition  the  interval 
Now  let 

n  n 

w\M)  =Y.Y.  WiiWMij  (12) 

»= 1  3- -1 


be  the  total  edge  weight  of  a  matching  at  time  tk  where  M. 
are  admissible  matchings.  The  total  edge  weight  wk(Ai) 
is  constant  on  the  interval  [tk,  ffc+i).  The  minimum  weight 
matching  is  calculated  at  the  beginning  of  each  time  interval 
and  remains  constant  over  the  interval  and  is  given  by 

n  n 

M(tk)  =  nmtwk(M)  =  min  EE  WjitkWij  (13) 

i= 1  3= 1 

The  total  weight  of  the  matching  over  the  interval  is  the 
weight  at  the  beginning  of  each  subinterval  time  the  length 
of  the  interval,  and  the  optimization  problem  becomes 

/-i 

Y]minu;fc(Af)(ffe+i -ffc)  (14) 

k—0 

f— 1  n  n 

=  E  “j*,11  EE  Wij(tk)Mij(tk+i  -  tk)  (15) 

k—0  i—1  j—1 

/-!_ 

=  ^M{tk)(tk+1-tk)  (16) 

k= 0 

(17) 


Next,  to  show  that  as  the  partition  gets  finer,  the  discrete¬ 
time  problem  approaches  the  continuous  time  problem,  we 
must  show  that  the  total  edge  weight  is  Riemann  integrable. 
Let 

n  n 

g(tk)  =  nun  E  E  ]V‘Afk)M,,  =  M(tk)  (18) 


i=l  j= 1 

/ 

I(g )  =  ^ErffcJfo.  “  tk~^ 

k—1 

(19) 

f 

=  lim  V  M(tk)(tk  ~  tfc-i) 
l!P|  |-°  f-r 

(20) 

where  ||P||  =  maxp  \tk  —  tk- 1|-  Since  A4(tk)  is  constant 
and  finite  on  the  interval  [tk,tk  +  1),  1(g)  exists  and  is 
finite.  Hence,  by  Riemann  integrability 

M(t)  =  [  mintn(Af)dt  (21) 

Jt i  M 

f 

=  hm  y^M(tk)(tk-tk-i)  (22) 

^P^_>0fe=i 

which  shows  that  in  the  limit,  the  discrete  time  case 
approaches  the  continuous  time  optimal  solution. 


IV.  Satellite  Application 

In  the  following  we  will  provide  additional  practical 
details  to  applying  the  Hungarian  method  to  the  problem  of 
interest.  Applying  this  algorithm  to  a  dynamic  system  in¬ 
volves  selecting  the  weight  functions  based  on  the  dynamics 
of  the  system  that  capture  the  behavior  to  be  optimized.  The 
dynamic  system  of  interest  in  this  paper  is  a  constellation 


of  low-earth  orbit  (LEO)  satellites.  The  function  of  these 
satellites  is  to  track  a  specific  location  on  the  earth  for 
communications,  imaging,  etc. 

We  consider,  the  problem  where  the  targets  are  not  nec¬ 
essarily  known  and  changing  with  time.  Another  dynamic 
aspect  of  this  problem,  that  arises  because  of  the  low-earth 
orbit,  is  that  targets  are  constantly  moving  in  and  out  of 
range  of  a  satellite  changing  the  edges  of  the  graph. 

We  start  by  defining  the  satellite  target  system  as  a 
weighted  bipartite  graph.  The  set  of  all  available  targets 
will  be  denoted  by  T(t)  and  the  set  of  all  satellites  will 
be  denoted  by  S(t).  Due  to  the  dynamics  of  the  problem 
each  satellite  will  only  have  a  subset  of  the  targets  in 
range  at  any  time.  Also,  at  any  given  time,  there  may  be 
satellites  with  no  targets  in  range  and  targets  that  are  not 
in  range  of  any  satellite.  These  targets  and  satellites  cannot 
be  matched  so  they  have  no  bearing  on  the  solution  and 
just  add  to  the  computational  load,  so  we  define  the  subset 
of  targets  T(t)  C  T(t)  as  the  instantaneous  set  of  targets 
{Tj,...  ,Tm}  that  are  in  range  of  a  satellite,  similarly  we 
define  S(t)  C  S(t)  as  the  instantaneous  set  of  satellites 
{Si, ,  Sm}  with  at  least  one  target  in  range. 


ST  ST  ST 


Fig.  5.  Removing  the  targets  and  satellites  with  no  edges  incident 


We  will  assume  that  there  are  i  =  1, . . .  ,n(t)  satellites 
and  j  =  1,.. .  ,m(t)  targets  in  S(t )  and  T(t)  at  time  t  £ 
[h,tf]. 

Next  we  want  to  describe  the  targets  in  T  that  are  in  range 
of  a  specific  satellite  in  S  at  any  time.  We  define  T!(t)  C  T 
as  the  set  of  targets  that  are  in  range  of  the  ith  satellite. 
Using  Fig  5  as  an  example  T1  =  {T)},  T2  =  {Tj,  72,13}, 
T3  =  {T1;T3},  andT4  =  {T2}. 

Note  that 


U  Tj(t)  =  T(t),  t  £  [t0,tf] 


(23) 


j  =  l,. ..,71 

The  definition  of  Tl(t)  allows  us  to  define  the  edge 
matrix 

'0  Tj  ^  T\t) 

1  Tj  £  T\t) 

where  i  =  1, . . . ,  n(t),  j  =  1, . . . ,  m(t).  For  Fig.  5  the  edge 
matrix  is 

10  0 
111 
10  1 
0  10 


E  = 


(24) 


The  edge  matrix  defines  the  set  of  valid  matchings  at  that 
instant  in  time.  We  will  define  the  weight  matrix  by 


Wij(t,E(t))  4 


fij(t ) 
oo 


Eij  =  1 

Eij  =  0 


(25) 


where  Wij  :  [ti,tf]  x  {0,1}  — >  R+.  For  this  satellite 
system,  we  want  to  maximize  the  time  that  the  satellites 
spend  looking  at  a  target  or  minimize  the  time  that  the 
satellite  spends  changing  its  attitude  to  track  a  target.  The 
time  that  a  satellite  takes  to  make  an  attitude  maneuver  is 
proportional  to  the  angle  that  the  satellite  must  slew  through 
which  we  will  denote  by  9(t). 

For  each  target  in  range  of  a  satellite  there  is  a  slew 
angle  to  look  at  that  target  based  on  the  target  location  and 
the  satellite’s  current  attitude.  We  will  denote  this  angle 
by  0ij(t).  Because  the  range  of  the  weight  assignment 
functions  is  the  positive  real  space  we  define  the  functions 
by 

fij(t)  =  {eij{t)f  (26) 

For  a  given  weight  matrix  the  total  weight  of  a  matching 
at  time  t  £  [f  i  ,  f /]  given  by 

n  m 

w{m)  =  y,YSPv(  *))2m«  (27) 

7=1  3= 1 

is  proportional  to  the  sum  of  the  total  angle  that  each 
satellite  must  slew  to  track  a  given  matching  of  targets. 
This  means  that  the  optimal  assignment  perfect  matching 
with  yield  the  matching  that  minimizes  the  sum  of  the  slew 
angles  that  the  satellites  must  move  through  to  track  those 
targets. 

Assuming  that  the  satellite  has  an  asymptotically  stable 
attitude  controller,  this  means  that  for  a  constant  matching 
the  total  matching  weight  is  going  asymptotically  to  zero. 
The  speed  with  which  the  satellite  can  acquire  a  target  can 
effect  the  algorithm  adversely.  If  the  attitude  control  system 
is  slow,  it  may  take  longer  to  slew  to  track  the  target,  then 
the  time  the  target  is  in  range.  Figure  6  shows  a  simulation 
of  how  the  total  weight  of  the  matching  varies  with  time. 


x  105 


time(sec) 


Fig.  6.  The  total  weight  of  the  matching  versus  time 


V.  Simulation  Setup 

In  the  following  we  will  discuss  the  simulation  of  the 
algorithm  on  various  LEO  satellite  constellations  used  for 


Fig.  7.  GlobalStar  satellite  consellation  screenshot  from  the  simulation 


sensing,  and  the  simplifying  assumptions  used  in  the  simu¬ 
lation.  The  simulation  was  used  to  validate  the  algorithm 
with  satellite  attitude  and  orbital  dynamics  added.  The 
simulation  was  also  used  to  look  for  computational  issues 
with  the  algorithm  and  test  situation  that  might  lead  to 
the  algorithm  becoming  “unstable”  for  example,  dithering 
between  two  targets  without  tracking  either.  Using  real¬ 
istic  parameters  and  orbits  is  given  less  emphasis  than 
providing  a  scenario  where  the  behavior  of  the  graph  is 
rich  and  complex  enough  to  thoroughly  test  the  algorithm. 
To  achieve  this,  some  of  the  test  orbits  were  chosen  to 
have  a  group  of  satellites  converging  on  a  region  of  the 
earth  with  multiple  targets.  Many  of  these  orbits  would 
be  impractical  to  launch  a  satellite  into,  but  all  orbits  are 
feasible.  The  algorithm  was  also  tested  on  various  global 
coverage  satellite  constellations.  The  orbital  propagation 
uses  standard  Keplarian  two  body  dynamics  and  accounts 
for  the  earth’s  oblateness  but  not  J2  effects. 

For  simplicity,  the  satellites  simulated  are  assumed  to  be 
axisymmetric  where  the  dynamics  are  given  by 


x(t)  =  I 


e 

</> 


+  u(t)  =  Ix(t) +u(t),  I  €  R3x3  (28) 


where  9,  and  V’  are  Euler  angles.  The  satellite  is  assumed 
to  have  a  maximum  slew  angle  giving  the  satellite  a  clearly 
defined  footprint  where  it  can  track  targets.  This  is  done 
to  prevent  the  satellite  from  looking  at  a  target  with  a  low 
angle  of  incidence. 

A  proportional  integral  attitude  controller  was  chosen 
with  gains  that  provide  an  average  slew  time  of  10  seconds. 
When  the  satellite  does  not  have  a  target  to  track,  it  tracks 
its  ground  track  to  return  the  attitude  to  the  center  of  the 
footprint. 

The  simulation  was  also  used  to  test  the  robustness  of 


the  algorithm  to  communication  failures.  The  algorithm  was 
implemented  in  a  distributed  fashion  which  improves  the  ro¬ 
bustness  of  the  algorithm.  In  the  event  of  a  communication 
failure,  the  behavior  of  the  algorithm  becomes  suboptimal, 
but  the  overall  objective  of  the  satellite  constellation  is  still 
accomplished.  In  all  cases  tested  the  algorithm  performed 
as  planned. 

VI.  Conclusion 

In  this  paper  we  investigate  the  optimal  assignment  of 
targets  to  LEO  satellites  with  narrow  beamwidths.  We  ex¬ 
tend  the  optimal  assignment  problem  to  dynamic  weighted 
bipartite  graphs  and  show  the  conditions  on  the  edge  weight 
assignment  functions  that  lead  to  the  optimal  solution  of  the 
assignment  problem  for  dynamic  graphs.  We  also  provide 
conditions  for  the  existence  of  a  perfect  matching  on  a 
given  bipartite  graph,  and  a  method  for  constructing  the 
smallest  virtual  graph  that  does  not  change  the  minimum 
matching  on  the  original  graph  and  guarantees  the  existence 
of  a  perfect  matching.  Finally  we  apply  the  algorithm  to  a 
constellation  of  LEO  satellites  and  verify  by  simulation  the 
effectiveness  of  the  algorithm. 

Future  extensions  to  this  work  include  investigating  the 
effects  of  time  delay  on  the  optimality  of  the  algorithm  and 
investigating  the  behavior  and  uses  of  different  edge  weight 
functions. 
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